The Effect of Curvaton Decay on the Primordial Power Spectrum 
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We study the effect of curvaton decay on the primordial power spectrum. Using analytical ap- 
proximations and also numerical calculations, we find that the power spectrum is enhanced during 
the radiation dominated era after the curvaton decay. The amplitude of the Bardeen potential is 
controlled by the fraction of the energy density in the curvaton at the time of curvaton decay. We 
show that the enhancement in the amplitude of the primordial curvature perturbation is, however, 
not large enough to lead to primordial black hole overproduction on scales which re-enter the horizon 
after the time of curvaton decay. 
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I. INTRODUCTION 



Recent results from experiments at the LHC in Geneva give a very strong indication of the existence of a Higgs like 
particle. If confirmed, this would be the first observation of a scalar field outside cosmology. For many years scalar 
fields have been nearly indispensable in theoretical cosmology, not only in solving many problems of the hot big bang 
model, but also in generating the primordial power spectrum of density fluctuations (see e.g. Ref. [l]). 
, The curvaton scenario [2|-y| is an elegant extension to the standard infiation scenario. Whereas in the standard 
Clj scenario the fields driving inflation and the field generating the primordial power spectrum are the same, in the 
Q , curvaton scenario these tasks are performed by two (or more) fields. While the infiaton field drives infiation, the 
• curvaton field is a mere spectator field (massless or nearly massless), picking up a nearly scale invariant spectrum during 
" inflation. After inflation the inflaton field decays into radiation and the curvaton behaves like dust, oscillating around 
. the minimum of its potential. After this brief matter dominated period the curvaton also decays into radiation. Due to 
its simplicity and elegance the curvaton scenario has enjoyed widespread attention in recent years (see e.g. Refs. ^5l-[l5|. 
I , In most standard curvaton calculations the power spectrum of the curvature perturbation Q (the quantity that 
J> ■ later on sources the CMB and LSS) is assumed to be directly inherited from the curvaton field. In this paper 
CNj , we study in detail the effect of curvaton decay on the primordial power spectrum. In particular we investigate 

■ whether the change in the initial power spectrum on small scales can lead to primordial black hole (PBH) pro- 
duction. To keep our results as general as possible, we do not specify a potential for the curvaton field, assuming 

^ ; ■ only that it oscillates around the minimum of its potential, and hence behaves like dust, before decaying into radiation. 

o\\ 

. The outline of the paper is as follows. In Sections HIl and Hill we present the equations governing the background and 

■ perturbation evolution. In Section IIVI we derive analytical approximations for the Bardeen potential, and calculate 

. . , the effect of curvaton decay on the power spectrum of C,. We then calculate the abundance of PBHs formed in Section 

^ ■ El We conclude with a brief discussion in the final section. 
• ^ ■ 

X: 

■ II. BACKGROUND DYNAMICS 

In this section we provide the background dynamics. We consider two interacting fluids in a Friedmann-Robertson- 
Walker (FRW) universe: radiation with equation of state P = uirPr, where Wr = 1/3 and the curvaton with dust like 
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equation of state. We assume the inflaton entirely decays into radiation at the end of inflation. Initially radiation 
dominates and the curvaton energy density is very sub-dominant. We assume that the curvaton decays into radiation 
with a constant decay rate T when the Hubble expansion rate H drops below F. 



A. Background Equations 

The continuity equations for each individual fluid in the background can be written as [l^ [13] 

Pa + SHpa = -Tpa- , (1) 

Pr + AHpr = Tpa , (2) 

where and pr are the background energy density of the curvaton and radiation respectively and a dot indicates 
the derivative with respect to cosmic time t. The right-hand side of equations ([T|) and ^ describe the background 
energy transfer per unit time from the curvaton and to the radiation fluid respectively. The Hubble parameter H is 
determined by the Friedmann constraint 

= jj^^ (P^ + Pr) , (3) 



where Mp\ = I/vSttG is the reduced Planck mass. In order to obtain the evolution of the above dynamical system, 
it is useful to write the equations in terms of dimensionless parameters 

n,^^, and a-—, (4) 



where 



^Pa^ P: (5) 



is the total energy density and hence ilcr + Sir = 1 • One can easily derive the following evolution equations for these 
dimensionless parameters: 

n',{N)-na{N)nr{N) + -^na{N) = o, (6) 

n',{N) + nr{N)na{N)-j^na{N) = o, (7) 

H'{N) + ^H{N) + ^H{N)nr{N)^0, (8) 

where a prime denotes differentiation with respect to the number of e-folding N = J Hdt. We set iV = at the end 
of reheating. 



B. Sudden decay limit 

Useful intuition can be gained by flrst considering the analytical solutions of the system of Eqs. ([6]-[8]) in the sudden 
decay limit where F = 0. This is a good approximation until the time of curvaton decay at = when T ^ H 
and the curvaton decays to radiation. In the sudden decay limit where T/H = 0, Eqs. ©, ([T]) and ([5]) can be solved 
easily for the background quantities ^r{N),i}cr{N), and H{N), to obtain 

N 

"^(^) - TT7^> (10) 

H{N) = i/oe-'~(l+pe^)'/' , (11) 
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where Hq is the value of H at the start of radiation dominated epoch, after the end of reheating, when N = 0. Here 
we have defined the ratio p 

P=-^ — > (12) 

where i„ and fir in are the initial values of VL^ and fir at the end of reheating, so that if radiation initially dominates 
p < 1. 

A key parameter in curvaton analysis is /d, the weighted fraction of curvaton energy density to the total energy 
density at the time of curvaton decay, defined by 

= +V ''^^^ " I^^ldcc - ^-Idcc • (13) 

As is well-known 0, Q, the non-Gaussianity parameter /nl on CMB scales is approximately given by /nl ~ l//d- In 
Section HVl we find the Bardeen potential 5' for various limiting values of /d. 



C. Beyond the sudden decay limit 



Here we provide solutions for ^li-{N) and ^lcr{N) which can also be used after the time of curvaton decay. Comparison 
with the full numerical solution to Eqs. ([SH5]) shows that while our analytical solution for H{N) in the sudden decay 
limit, Eq. (jlip is in reasonable agreement with the numerical solution the analytical sudden decay solutions for VL^{N) 
and ^a{N), Eqs. © and (ITUl) . are not accurate for N > N^as expected. This indicates that more accurate analytical 
solutions for Vt^{N) and Vla{N) can be found by inserting the solution for H{N) into Eqs. ^ and ([7]) and then solving 
for i},{N) and rJ^(A^). 

Defining 



from Eqs. ^ and ([7]) we find 



(14) 



X' = X--{l + X)X, 



(15) 



with H{N) given by Eq. (fTTI) . It is not possible to solve Eq. ([T5|) in general. However, away from the time of curvaton 
decay, we have X < 1, i.e. fl^ < so Eq. (jlSp can be approximated by a linear differential equation which can be 
solved analytically yielding 



X 



pexp 



pexp 



N 



dN' 1 



2T 



H{N'] 



pe 



N 



l+pe 



N 



1 



N 



(16) 
(17) 



In the integrand above H{N') is calculated from the sudden decay limit solution given in Eq. (jlip . Correspondingly 

1 „ X 



l + X 



i + x 



(18) 



Note that before the time of curvaton decay when one can neglect the last term in Eq. ([15]) and X ~ pe , Eq. 
reduces to Eqs. (O and ([TU|) obtained in the sudden decay limit. 

For the physically interesting case in which the curvaton makes a sub-dominant contribution to the total energy 
density at the time of decay, corresponding to f2cr|doc ^ 1, one can take ilcr{N) <C i^i{N) throughout the whole 
evolution and to a good approximation H{N) is given, as in a radiation dominated Universe, by H{N) ~ Hqc^'^ . 
Inserting this expression for H into Eq. (|16p yields 



'- p exp ( X 



r 



-IN 



(/d « 1) 



(19) 
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FIG. 1: The evolution of fir and fla for F/Hq — 10~*. The left panel is for /d, the weighted ratio of the curvaton energy 
density to the total energy density at the time of curvaton decay, equal to 0.08 and the right panel for /d = 0.16. In each plot, 
the upper curves are Qr while the lower curves are Q^. The solid red curves are the full numerical solution, the blue dotted 
curves the analytic solutions, Eqs. (jlSp . and the green dashed curves the analytic solutions in the /d ^ 1 limit, Eqs. H19p . As 
expected, for small values of /d the agreement between the numerical and the analytical solutions is very good. 

In Fig.[T]wc have plotted the full numerical solutions for ^^{N) and ^1,-{N) and compared them with our analytical 
solutions Eqs. and As can be seen, for small values of rialdcc, or equivalently small values of /d, the 

agreement between our analytical solutions and the full numerical solutions is very good. 

Finally we obtain an estimate for N^, the time of curvaton decay. A good criteria to define in the sudden decay 
limit is when the last term in Eq. becomes comparable to the second term and X'{Nd) = 0, and the ratio flcr/flr 
reaches its maximum. This gives r(l + X{Nd)) = H{Nd)- In the limit where X{Nd) ^ 1, or equivalently /d <C 1, 
this reduces to the standard result that H{Nii) ~ F or 



(/d « 1) 



Inserting this expression for A^d into the definition of /d in Eq. (|T^ yields 

1/2 



, 3p (Ho 



(/d « 1) 



(20) 



(21) 



In general, when /d ^ 1/2 then Eq. (|20)) receives corrections and one has to find by solving X' — with X given 
by Eq. (HID. 



III. PERTURBATIONS 



In this section we study the perturbed Einstein and fiuid equations. The perturbed metric line element is |18l - l20l | 

ds^ = -(1 + 2<j))dt^ + 2aB jdMx' + [(1 - 2?A)% + 2E^,j] dx'dx' . (22) 
The perturbed Einstein equations are then 



and 



vI'' + ^* + ^* + _(i+.)C = 0, 



(23) 



(24) 
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where the time-dependent equation of state w is given by 



w . 



P 
P 



E 



Pa 
P 



(25) 



in which Wq is the equation of state for each fluid given by uja — Pa/ Pa with lo^ = and ujr = 1/3. 

Here we have defined the Bardeen potential, or curvature perturbation on uniform shear hypersurfaces, as 

* = V - H{B - E) , (26) 

and the curvature perturbations on uniform density slices C, and on comoving hypersurfaces, TZ, respectively as 

H6p 



p 



n = %i>-HV , 



(27) 



where V = a{v + B), and v is the total scalar velocity potential. 
Equations ([^5]) and ([M]) can be combined to give 



(28) 



In particular, we see that on large scale where k/aH 0, C — —T^- 
The equations of motion for each fluid are 

P 1 

Sp'^ + 3{5pa + Spa) - 3{pa + PaW TTf^Pa + Po^iVa + CF s) - -;^{Qa4> + ^Qa) 



V 



Qa 



H{pa +Pa) 



(1 + cl) - 34 



1 



H H{p^+pa,) 



[Spa-QaV] = 0, 



(29) 
(30) 



where the sound speed for each fluid is defined by = Paj Pa- We shall assume further below that each fluid is 
intrinsically adiabatic so Pa = Pa{pa) and 6pa = c^Spa- Also Va is the scalar velocity potential for each fluid and 
14 =a(wa +B). 

In order to express the fluid equations, Eqs. (j29D and pop. in gauge invariant form, we define the curvature 
perturbations (a and TZa for each fluid as 



Ca = -W- H^- 
Pa 



and 



One can cast the perturbed fluid equations, Eqs. (j29[) and into [iJ 



Ca-- 



H'Qa / 5pa_ _ Sp 
H'p'a V P'a P 



7)-Jh(sQ<'-^sp^+ 

' J Hp'^ V Pa J 



^ - [1 



Wa 



and 



Qo 



H' 



Hipa+Pa) H 



(7^-7^„)- 



Pa 



Pa +Pa 



cl (7^„+Ca) . 



(31) 



(32) 



(33) 



(34) 



In the absence of energy transfer between the fluids, Qa — SQa — 0, we see that (a for each fluid is constant on 
super- horizon scales [2l[. Note that in deriving Eqs. ([55)) and we have assumed that each fluid is intrinsically 
adiabatic so pa = Pa{Pa) and 5pa = c^Spa- 

We also note that C and TZ defined in Eq. (|27|) can be written as the weighted sum of (a and TZa, 



(I +UJa)Pa 

+ uj)p 



Ca J 



(35) 
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and 



So far we have not specified the perturbations in the energy transfer hQa- FoUowing Ref. [1^1 , we assume that the 
decay rate F is fixed by the microphysics so 5F = and therefore 

^Qa = -^Spa , 5Q^ = Tdpa ■ (37) 

The system of equations in terms of {Ca,T^a, Ci ^} is given by 

C = ;^77^^(C-C.) + :T-7TH*--7IHh7T^P^- (38) 



2(F + 3i/)' ' ia'^m \3H + T 

F , H'\ 3(1 + f7^). 



Note that this is a closed system of equation for {CcTi-a, CjT^} that ^ can be ehminatcd from these equations 
using Eq. in terms of C and TZ. 

Alternatively, one may write the system of equations in terms of {Cr, T^-r, C: ^} or {C-cT^r, Ca,TZa} as given in 
Appendix El 



IV. ANALYTIC CALCULATION OF BARDEEN POTENTIAL 



In this section we provide the analytical solutions for the Bardeen potential ^ in different limits. As can be seen, the 
system of Eqs. ([55)) -(PT |) is too complicated to be handled analytically for all modes. However analytical solutions for 
'I' can be obtained in some limiting situations. In the next two subsections we consider modes which are super-horizon 
at the time of curvaton decay, k < a{Nd)H{Nd), and modes which are always sub-horizon during curvaton evolution 
corresponding to k > Hq where Hq is the Hubble constant at the end of reheating when a{N = 0) = 1. 



A. Super-horizon modes 

Here we provide the solution for the modes which are super-horizon at the time of curvaton decay and re-enter the 
horizon during the second radiation stage. First we consider the epoch before the curvaton decays, N < N^- In the 
sudden decay limit, for the super-horizon modes from Eqs. (jAl[) and (jA2p it can be shown that Cr — and ~ 0. 
As expected, on super-horizon scales both and Cr remain frozen so one can approximate ^o- with its value at the 
time of horizon crossing during inflation 

UN) ^ U . (42) 

Also to further simplify the analysis, we consider the conventional curvaton mechanism where C^in = and there is 
no initial radiation perturbation, corresponding to entropic initial conditions. In this limit, either by solving Eq. (jA7p 
or using the definition (j35p . we have 

Substituting this into Eq. ([25)1 and noting that w = ili-{N)/3 and using Eq. PT|) . results in 



2 ^ 2(l-fpe^) ^ ' 



7 



This can easily be solved with the result 



C y/l+pe^ _ ^ (16 + 8p - e 



-3N 



(45) 



Here C is a constant of integration which is obtained by matching this solution to the value of ^f" at the end of inflation, 
which gives 



C 



1 



*(0) + 



3C^ 



16 + 8p - 2p^ + 



5p'^ 



(46) 



To obtain the second approximate relation we used p <C 1 and ^'(0) ~ which is a good approximation for all modes 
at the end of inflation [22| • 

Having obtained ^' during curvaton evolution, we now find after the curvaton decays. The governing equation 
for 4" during the radiation domination stage with Wr — 1/3 has the standard form 



fc2 



= 0. 



3a2iJ2 
Using 

during the radiation era this leads to a solution in terms of Bessel functions, 

^(N) = [ci J3/2 (fce^) + ^2^3/2 (ke^)] , 

in which ci and C2 are two constants of integration and 

_ fce-2jVa 

" V3i/d ■ 



(47) 

(48) 
(49) 

(50) 



In the limit where the curvaton makes a sub-dominant contribution to the total energy density at the time of decay, 
corresponding to /d <C 1, one can use Eqs. (ITlT) and ((20)) to obtain — F and k ~ k/Hg. Also the condition for the 
mode k to be superhorizon at the time of curvaton decay, k > Hde^"^ , translates into k > ^/THq. 

We can now determine the constants of integration ci and C2, by requiring that both and VP' are continuous at 
the time of curvaton decay, N ~ Nd- This gives 



Cl 



C2 = 



fce^^^/'n'/2(^e^^) - |e3^^/2y3/2(fce^^) W(iVd) - e^'^'^/^Y./^ike''-)^' {N^) 



-ke^^-'^4,^{ke''-) + ^e^^^/V3/2(fce^^) *(iVd) + 6^^-/^3/2 (fce^^)*'(iVd) 



(51) 
(52) 



Note that in the above expressions ^'(A^d) and 'I''(A^d) are calculated from the solution obtained from the period 
before the curvaton decay, Eq. dJS]). 

By construction, we know that a{Nd)H{Nd) < k, so ke^"^ ~ ke"^"^ / < 1. Using the small argument limits of 
the Bessel functions we find that C2 <C ci and therefore 



vI/(iV)~cie-37V/V3/2(fce^), iN>Nd). 



(53) 



Now suppose these modes (with ke < 1) re-enter the horizon during the second radiation era at the time N — Nk- 
The value of can be estimated by calculating when the argument of the Bessel function in Eq. (15^ becomes of 
order unity and Jz/2{ke^) starts to oscillate. This yields ke^'' = 1 or Nk = — Ink. During the period Nd < N < N^ 
in which the mode is still super-horizon during the second radiation stage we can use the small argument limit of the 
Bessel function Ji,{x) ~ {x/2Y /T{iy + 1) (here T{i/ + 1) is the gamma-function) to obtain 



*(7V) 



Cl 



r(5/2) V2 



3/2 



Nd<N <Nk 



(54) 



This is a very interesting result; for A^d < N < Nk, 'J is constant with the value given by Eq. ([5^ . Numerical 
evolution of 4*, shown in Fig. [2j verifies the existence of this plateau. 
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FIG. 2: The evolution of the Bardeen potential ^ for super-horizon modes with fd ~ 0.08 and 0.16 (left and right plots 
respectively). In both plots F/Ho — 10~* and the red curve is the full numerical result while the black dotted curve is our 
analytical solution, Eq. (|49[) . In both plots, the existence of the plateau, Eq. (I54|) . for the modes in the range A*'d < N < Nk is 
evident. As expected from Eq. H58p . the larger /d is, the larger the amplitude of ^ is. 



During the epoch N > Nk, after the mode k has re-entered the horizon, the Bessel function in Eq. (|53)) oscillates 
rapidly. Using the large argument limit of the Bessel function we obtain 

*(iV) ~ ,/Xe-2JV^^ cog (^^N > (55) 

V TTK V 4/ 

As we shall see in the following section in order to study whether PBHs are overproduced, we need to estimate 
^'(A'd) and ci. As mentioned before, the key parameter in controlling the amplitude of ^'(iV) during the second 
radiation era is fd- Here we calculate ci for two extreme cases (i): /d ^ 1 corresponding to pe^'^ <C 1 and (ii): 
/d > 1/2 corresponding to pe^'^ <C 1. 

For the limit /d ^ 1, from Eq. (|45p we obtain 

vI/(iVd)=^-^pe^^C.. =^-yC.., (/d«l), (56) 

where to obtain the second approximation, Eqs. (|20p and (piT) have been used. Similarly, using Eq. to eliminate 
^' we obtain 

*'(iVd)~-^C.. ~vi/(iVd), (/d«l). (57) 



As a result 



ci ~ |fc"3/2 ( 3^(^^) _^ ) ^ - V2^/dfc-3/2c,. , (/d « 1) . (58) 



We have checked this numerically and our analytical estimation of ci is in good agreement with its numerical value. 
On the other hand, in the limit where pe^'^ 3> 1, so fd ^ 1/2, one finds 

vI/(iV,)^-^C. , vi/'(iV,)^_^vi/(iV,)-^C.^O. (59) 

Inserting these results into the expression for ci in Eq. ([5T|) yields 

This indicates that for large enough /d, the amplitude of during the second radiation dominated epoch is nearly 
independent of fd- 



FIG. 3: The evolution of the Bardeen potential ^ for sub-horizon modes with k/H{0) = 2 and 10 in the left and right panels 
respectively, for /d — 0.16 and F/Hq = 10~*. In both plots, the thick dashed blue curve is the full numerical result while the 
thin solid red curve is the numerical solution obtained from Eq. (|61|l . 



In Fig. [2] we have plotted ^'(iV) for the super-horizon modes with /d ^ 1 and /d > 1/2. The main result here is 
that the amplitude of 5* increases during the final radiation dominated era due to curvaton dynamics. However, the 
increase in amplitude of ^ is less efficient for small /d. We shall see whether this will have implications for primordial 
black hole formation in Section IVl below. 



B. Sub-horizon modes 



In this sub-section we calculate 5* for the sub-horizon modes. In general, it is not easy to solve the system of equation 
in this limit. A particular limit which may be handled semi-analytically is the case where the term (k/aH)'^ ^ 1 
throughout curvaton dynamics, corresponding to fc > Hq. These are modes which are sub- horizon during the entire 
curvaton dynamics. 

Differentiating Eq. and using Eq. ([5^ and Eq. (IMl) to to eliminate C' and TZ we obtain the following second 
order differential equation for 4* 

2{A-n„) 3(4 - r!,) a2i72 " ^" ^ 

In the limit — 0, this reduces to the standard equation for in a radiation dominated background given in 
Eq. gZD. 

Equation (pT|) can not be solved analytically, as far as we know. In Fig. [3] we have plotted 'i> for different values 
of k corresponding to modes which are deep inside the horizon at the end of inflation. In principle one could solve 
Eq. (|5T|) semi-analytically and find the values of vE'(A^d) and ^'{N^) and insert them into Eq. ([55]) to find \1/ in the final 
radiation dominated era. With an analytical solution for in this regime one could then calculate the abundance of 
PBHs formed from sub- horizon fluctuations, c.f. Refs. [13: H^l- This is beyond the scope of this work however, and 
we focus in the next section on the 'standard' case of PBHs forming when super-horizon modes reenter the horizon, 
using the analytic solution for ^E* from subsection IIV Al in the final radiation dominated era, i.e. Eq. ([5^ , with ci 
given by Eqs. dSH]) and pi)) . 
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V. PRIMORDIAL BLACK HOLE FORMATION 



Primordial black holes (PBHs) are a powerful tool for constraining models of the early Universe. Due to their 
gravitational effects and the consequences of their evaporation there are tight constraints on the number of PBHs that 
form (see e.g. Refs. (H,!!^)- These abundance constraints can be used to constrain the primordial power spectrum, 
and hence models of the early Universe, on scales far smaller than those probed by cosmological observations. We 
found in Sec. IIVI that the power spectrum is enhanced during the radiation dominated period after curvaton decay. 
In this section we therefore investigate whether this enhancement is sufficiently large to lead to PBH over production. 

A region will collapse to form a PBH, with mass of order the horizon mass at that epoch, if the smoothed density 
contrast, in the comoving gauge, at horizon crossing {R = Cs[aH)^^ where Cs — l/VS is the sound speed), ScomiR — 
1/aH) = Shor{R), satisfies the condition Shor{R) > 5c [11], where 6c ^ 1/3 ^. The fraction of the energy density of 
the Universe contained in regions dense enough to form PBHs is then given, as in Press-Schechter theory [28j . by 

/? = r P(4or(i?)) d4or(i?) . (62) 

Assuming that the initial primordial perturbations are Gaussian, the probability distribution of the smoothed density 
contrast, P((5hor(-R)), is given by (e.g. Ref. [J) 

1 ..^r «1o,(R) 



^<*"<«»^7ia'--'l-^j- '''' 

where a{R) is the mass variance 

cT^{R)= / W^ikR)rs{k,t)^. (64) 

Here W{kR) is the Fourier transform of the window function used to smooth the density contrast, which we take to 
be Gaussian, so that W{kR) — exp (— A:^i?^/2) and Vs{k,t) is the power spectrum of the comoving density contrast, 

Vs{k,t) = ^{\Scon.{k,t)\') . (65) 

Inserting the expression for the probability distribution, Eq. (|63p. into the Press-Schechter expression for the initial 
PBH abundance, Eq. (|62|) . gives 

/3 = ^ / exp f-T^MzT^l d(5hor(i?) = erfc ( ^ . (66) 

V2^ahor(i?)Ae V 2aUR)J '^"'^ ^ \V2ai,or{R)J ^ ' 

The constraints on the initial PBH abundance are translated into constraints on the mass variance by inverting this 
expression. The constraints are scale dependent and lie in the range 13 < 10~^° — 10~^ [U, The resulting 

constraints on crhor(-R) lie in the range ahoi-{R)/Sc < 0.1 — 0.2 [23 |. 

In order to calculate (Thor(P) we need to evaluate the density contrast at the epoch when the scale of interest, 
R^^ = kc, enters the horizon, Cgfcc = aH, at TV = Nk- We consider PBH formation for modes which re-enter the 
horizon after curvaton decay, kc > a{Nd)H{Nd), for which we can use the analytical results from section HV Al 

Using the Poisson equation and Eq. I^E^ the comoving density contrast is related to the Bardeen potential by 

As a result the number of PBHs formed is controlled by the amplitude of ^; the larger the amplitude of ^, the larger 
the abundance of PBH formed. As we saw before, the amplitude of 4* decreases as /d is decreased. Let us consider 



^ It was previously thought that there was an upper limit on the size of fluctuations which form PBHs, with larger fluctuations forming 
a separate closed universe, however Kopp et al. [27II have recently shown that this is in fact not the case. 
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the case where the PBH formation is most efficient, corresponding to fd ^ 1/2. Using Eq. (p7|) and Eq. ([55)1 for the 
Bardeen potential evaluated at iV = Nk with the constant ci given by Eq. 

- 4 /_^^' 



^^C.f-^V(^e-^)-^[^3/2(^e-^r • (68) 



25 Va^f, 

Using Eq. (|48)) the time of horizon crossing, Nk, is related to the time of curvaton decay, N^, by 



so that k, defined in Eq. ([50)) . is given by 



ke""' = ^ . (70) 



Inserting the expression for the power spectrum of the comoving density contrast, Eq. (|68p . into the definition of the 
mass variance, Eq. ([M]) . and using Eq. ([70]) gives 



25c4 

Finally, using the numerical approximation 



r^e-'^'/'' [J3/2(fc/fcc)] 

Jo "'c 



^loAR) = ^-Pcl :pe-'=-/'=^ \J3/2{k/k,)r . (71) 



dxe-^ J3/2(a;)^ Ci; 0.02, (72) 

"'0 

we find 

al,{R)c,OAVc^. (73) 

Therefore the requirement to avoid PBH overproduction, o'hor(^)/'5c < 0.1 — 0.2, leads to a straight-forward, and 
fairly weak, constraint 'P(^^{R) < 10^^ — 10^^ which is easily satisfied for an almost scale invariant curvaton field. 

On the other hand, if we consider the limit in which /d <C 1 the above result becomes Choi(^) ~ fd'^Ca- ^ 
result, the condition on Vq^ [R) becomes even less restrictive as expected. 



VI. CONCLUSION AND DISCUSSION 



In this paper we have studied what effect the curvaton decay has on the primordial power spectrum of the density 
fluctuations and the Bardeen potential. To this end we studied a simple system comprising only radiation and the 
curvaton, which we modelled as a pressureless fluid, using a flat FRW universe as background. 

The key parameter in our analysis is the weighted fraction of the curvaton energy density to the total energy density 
at the time of curvaton decay, /d, defined in Eq. p^ . Solving the system of governing equations analytically, we 
found that on super-horizon scales an increase in /d will lead to an enhancement of the amplitude of the Bardeen 
potential due to curvaton dynamics. Unfortunately we were not able to solve the system in the small scale limit 
analytically, and therefore leave semi-analytical solutions in this regime to future work. 

Having established an enhancement in the density contrast on super-horizon scales, it is natural to ask whether 
this increase will lead to observational consequences, in particular to the overproduction of PBHs. We studied this 
issue in detail and found that the enhancement is too small to lead to significant PBH production. We can therefore 
conclude that the enhancement of the primordial power spectrum on super-horizon scales does not lead to additional 
constraints on the curvaton model through PBH production. 

However, since we only found analytical solution for the super-horizon modes, we could not investigate whether 
there would be significant PBH production on sub- horizon scales. As shown in Refs. [13] and [1^, PBH production on 
sub-horizon scales can have an significant effect leading to further constraints on the model. We hope to investigate 
these questions in future work. 
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Appendix A: Fluid equations in different forms 



The closed system of fluid equations in terms of variables {C, 72., Co-, T^a- } is given in Eqs. 
the equivalent systems of equations in terms of {CR-XnT^i } and {Cr,T^iX<T,T^a }■ 
In terms of {C,7?., CtiT^i } the system of equations read 





r Tn„ 




6H r 


c = 


r 


K = 




K - 





-s 



T + 3H 

{3H + r){n, 



Vila - 4ff a 



Is- 



H 



a^m \3H + T 



Tin 



IT 
ir_ 



An, 



3n„ 



3 + a- 



3 + a- 

Here the entropy perturbation 5* is defined via 



3 + a 



5 = 3 - C ) = '3H 



Spa 



Pr 



(|41l) . Here we present 



(Al) 
(A2) 
(A3) 
(A4) 



(n, + Cr) 



(A5) 



We see that Eqs. (|Aip - (|A4p give a closed system of equations for four variables {Ca,T^a}- Note that vj/ can be 
eliminated in these equations from Eq. d^B]) in terms of ( and TZ which are expressed in terms of Cq and TZa from 
Eqs. ([55)) and ([55]) . As mentioned in Ref. there is an apparent singularity in the system above when p,- = and 
4Hil,. In order to overcome this problem it is convenient to trade Cr and TZ, for C and TZ as we did in Eqs. 



381) - (|4T|) in the main text. 
Alternatively, the system of equations in terms of {Cr, T^r, C, 7^} is written as 



c = 



r + 



TZr; 



TZ', = 
TZ' = 



H{3 + n,) 
fSTQ^ H' 

H' AHVLr 



H H{3 + n,) 



{TZ - TZ,) 
TZ 



1 - 



(TZ, + Cr) 



1 , ^ , , AHil, 

2(^ + "-)^ + ^/(fTT^ 



■Cr 



3a27J2 



(A6) 
(A7) 
(A8) 
(A9) 
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